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ABSTRACT 

Using a sample of galaxies from the Sloan Digital Sky Survey spectroscopic 
catalog with measured star-formation rates (SFRs) and ultraviolet (UV) pho- 
tometry from the GALEX Medium Imaging Survey, we derived empirical linear 
correlations between the SFR to UV luminosity ratio and the UV— optical colors 
of blue sequence galaxies. The relations provide a simple prescription to correct 
UV data for dust attenuation that best reconciles the SFRs derived from UV and 
emission line data. The method breaks down for the red sequence population as 
well as for very blue galaxies such as the local "supercompact" UV luminous 
galaxies and the majority of high redshift Lyman Break Galaxies which form a 
low attenuation sequence of their own. 
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Introduction 



Although directly tracing young massive stars, the ultraviolet (UV; A ~ 912 - 3000 A) 
luminosity of a galaxy is not a straightforward measure of its current star-formation rate 
(SFR), nor is in fact any other observable related to new born stars. In particular, the 
galaxy's dust content and past star-formation history (SFH) have a significant influence on 
the interpretation of the observed UV flux in terms of current star production. A large, 
sometimes dominant fraction of the UV emission may be obscured by dust and reprocessed 
at far inf rared (FIR) wavelen gths. The UV spectral slope is commor ily used to estima te this 
fraction (IMeurer et al.lll999l ) but it is strongly affected by the SFH (IKong et al.l 120041 ). Also 
because UV emitting stars live long enough for successive generations to coexist, the SFH 
over the past few hundred Myrs must be known to translate the dust corrected UV flux into 
a more instantaneous SFR, such as derived from the galaxy's Ha emission. In the case of a 
constant SFR and in the absence of dust, the UV luminosity to SFR ratio reaches a plateau 
after ~ 10® yrs (e.g. Kennicutt 1998), but a strong starburst will cause the UV luminosity 
to scale differently with the SFR. Interpreting the UV em ission of early-ty pe galaxies is also 
less straightforward due to contamination by older stars ( IRee et al.l 120071 ). 



Dust obscuration and SFH may be estimated with the help of additional data (e.g. the 
Balmer decrement or far- infrared emission for the dust, the Balmer break for the SFH) and 
of theoretical assumptions. However such additional information is not always available or 
in fact acquirable, in particular at high redshifts where SFRs are generally derived from UV 
and/or infrared (IR) photometry with rather large uncertainties. New near-IR spectrographs 
on 8-lOm telesco pes are now making it possible to detect Ha emission and/or continuum 
breaks at z > 2 Jvan Dokkum et al.l [2005l : lErb et all BoOGbl : iKriek et all boOGal lbl) but the 
technique is still limited to p rominent features . Detecting Ha at z ~ 2 implies a SFR 



greater tha n a few MfD vr ^ (Erb et al.l l2006bl ). which is not representative of the whole 



population (IKriek et al.ll2006bl ). Optical images which pick up the rest-frame UV at z > 2 
remain the easiest data to obtain. 

Here we use medium deep UV photometry from the Galaxy Evolution Explorer (GALEX) 
and the wealth of additional data provided by the Sloan Digital Sky Survey (SDSS) to de- 
rive simple empirical relations between the observed UV luminosity and the SFR of local 
star-forming galaxies. SFRs were derived for tens of thousands of SDSS galaxies using their 
emission lines and state-of-the-art models including a consistent treatment of the dust from 
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the UV to the far-IR (iBrinchmann et aLl 12004 ). We assume these SFRs to be the best pos- 
sible estimates at the present time, given the quahty of the spectro-photometric data and 
the tested rehabihty of the models, and present an empirical method to recover them from 
the UV luminosity of galaxies using their UV— optical colors. We compare our relations 
with existing methods, in particular attenuation estimates based on the slope of the UV 
continuum which are commonly used, and investigate their limitations both locally and at 
high redshift. 

The data are summarized in Section 2. In Section 3 we review the relation between the 
UV luminosity and the SFR of star-forming galaxies as well as several published methods for 
estimating their UV attenuation. In Section 4 we present simple empirical color relations that 
best reconcile the UV data with the SFR estimates based on emission line measurements, 
and discuss their limitations. Our conclusions are presented in Section 5. Throughout the 
paper we assumed a flat AC DM cosrnology with Hq = 70 km s^^Mpc"^, Qm = 0.3 and 
= 0.7, and a Kroupa IMF jKroupalboOlh . 



2. Data and derived physical quantities 



We select galaxies from the Sloan Digital Sky Survey (SDSS) spectroscopic catalog (Data 
Release 4; Adelman-McCarthy 2006) with NUV and FUV photometry from the GALEX 
Medium Imaging Survey (Internal Release 1.1; Martin et al. 2005, Morissey et al. 2005, 
2007). The UV filters have effective wavelengths of 1528 and 2271 A respectively. The 
Medium Imaging Survey (MIS) has a 5 cr detection limit of 22.7 (AB magnitude) in both 
filters for a typical exposure. This magnitude limit corresponds to a cut in magnitude error 
of ~ 0.1 in the NUV band and ~ 0.2 in the FUV band (Bianchi et al. 2007, their Fig. 2). 
Our primary sample consists of 23400 SDSS galaxies with r < 17.8, z > 0.005, measured 
Ha emission, aperture corrections less than 1.3 dex (defined as the ratio of the total SFR 
to the SFR estimated within the fiber; see below), and GALEX coverage in the NUV band. 
Adding FUV coverage reduces the sample to 17500 galaxies due to occasional failures of the 
GALEX FUV detector. Galaxies flagged as AGNs in the SDSS MPA/JHU DR4 value-added 
catalog^ have been excluded. 



The p hysical properties of S DSS galaxies were analyzed i n detail by iKauffmann et al. 



(l2003al lb[):[ Tremonti et al.l (120041 ) and IBrinchmann et al.l (120041 ) (hereafter B04) among oth- 
ers. In particular, the full likelihood distributio ns of their SFRs were deriv ed by fitting 
all strong emission lines simultaneously using the ICharlot &: Longhettil (l200ll ) models, fol- 



^http: / /www. mpa-garching.mpg.de/SDSS/ 
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lowing the me t hodol ogy of ICharlot et al.l (120021 ) (B04). Dust is accounted for with the 
Chariot fc Falll (120001 ) multicomponent model which provides a consistent treatment of the 
attenuation of both continuum and emission line photons. The dust attenuation is based on 
the Ha/H/3 ratio to first order but is really constrained by all the lines. B04 also devised a 
method for estimating the SFR of early-type galaxies with no detectable Ha emission from 
their 4000A break index but we excluded those from our sample. We use the medians of the 
SFR distributions and consider these values, noted SFR^ (for emission lines following B04), 
to be the best currently available estimates of the SFR given the quality of the data and the 
technique used to derived them. Uncertainties are discussed in detail in the original paper. 



Other quantities such as 4000 A break indices and stellar masses ( Kauffmann et al. 



2003al ) are also available from the SDSS MPA/JHU DR4 value added catalogs^. The stellar 
mass is defined as the total mass of stars formed over the lifetime of the galaxy. The 4000A 
break index - D„(4000) - is defined as the ratio of the average flux d ensity in the narrow 
bands 3850-3950A and 4000-4100A following [j^logh et all JlQQsh . It is a relatively dust 
insensitive measure of a galaxy's SFH, equivalent to th e ratio of the SFR av eraged over the 
last ~ 10^'^ years to the SFR averaged over > 10^ years (1 Johnson et al.ll2007l ). The D„(4000) 
distribution is strongly bimodal around D„(4000)~ 1.6, dividing galaxies into the well known 
'red sequence' of early-type, old star dominated galaxies (Dn (4000) > 1.6) and the 'blue cloud ' 
of late-type galaxies with recent star formation (D„(4000)< 1.6) (iKauffmann et al.ll2003al : 
Strateva et al.ll200ll ). In the following we refer to these 2 populations as simply red and blue 
galaxies. 

Our final 2 samples consist of the 20800 galaxies in the primary sample (89%) that have 
been detected by GALEX in the NUV band (SDSS-I-NUV sample), and of the 14900 galaxies 
in the primary sample with additional FUV coverage (85%) that have been detected both in 
the NUV and FUV bands (SDSS+NUV+FUV sample). In the following we use the larger, 
SDSS+NUV only sample whenever FUV fluxes are not explicitly needed. These samples 
are strongly biased against red sequence galaxies but complete for blue galaxies: ~ 98% of 
D„(4000)< 1.6 galaxies are detected in the NUV band (and in the FUV band when both are 
available) against ~ 72% of D„(4000)> 1.6 galaxies (~ 52% in both UV bands when both 
are available). The average magnitude error is 0.03 in the NUV band and 0.07 in the FUV 
band for the blue population; 0.08 in the NUV band and 0.16 in the FUV band for the red 
population. 

We derive absolute magnitudes in all the bands from the redshift and the Galac- 
tic extinction-corr e cted SDSS-I-GALEX photometry using the kcorrect v4_l software of 
Blanton fc Roweij (120071 ). In order to minimize the uncertainties on the fc-corrections, the 
magnitudes are /c-corrected to the mean redshift of the SDSS sample (z = 0.1) and are 
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noted ^■^mag where mag = f or n for the GALEX FUV and NUV bands, g, r, i and z 
for the SDSS bands. The fc-correction at redshift 1 is by definition —2.5 log(l + z) in all 
bands for all galaxies and deviates from this value towards both ends of the redshift range 
(0.005 < z < 0.28). In the UV bands, this deviation is less than 0.1 magnitude for 95% of 
the galaxies. 

Figure [1] shows the distribution of the SDSS+NUV sample in the "'^(n — r) vs D„(4000) 
plane. The '^■^{n — r) color distribution is strongly bimodal (Wyder et al. 2007) with 
0.1 — r) = 4 defining roughly the same boundaries as D„(4000)= 1.6 between the red 
and blue populations. Galaxies with ^-^{n — r) < 4 and D„(4000)< 1.6 (the 'blue cloud') 
represent 70% of the sample (81% of the SDSS+NUV+FUV sample); galaxies with "-^(ra — 
r) > 4 and D„(4000)> 1.6 (the 'red sequence') represent 22% of the sample (10% of the 
SDSS+NUV+ FUV sample). The so lid lines are polynomial fits to the dust/color/SFH rela- 



tion derived by I Johnson et al.l (120071 ) (see Section 3) for given values of the FUV attenuation 
as marked in the figure. The fits are good for galaxies with D„(4000)< 1.6 but less reliable 
for red sequence galaxies (see Johnson et al. 2007 for details). The model illustrates how 
broad band colors depend on both the SFH and the amount of dust attenuation. 



3. Deriving a Star Formation Rate from an Ultraviolet flux 

3.1. Calibration 



The SFR measured from the UV emission is usually written as: 



where r]uy converges to rfljy for a constant SFR. Scaled to a Kroupa IMF, the most commonly 
used factor is logirj'fjy) = 28.02 flKennicuttlll998r ). It assum es that the UV spectrum is nearly 



flat in L,y over the wavelength range 1500-2800 A. Using the lBruzual fc CharlotI (120031 ) stellar 
population synthesis models with similar assumptions (solar metallicity, a constant SFH and 
a Kroupa IMF), S. Salim (private communication) derived slightly higher factors for the 
GALEX filters: log{7]%y) = 28.09 and \og{7]%uv) = 28.08. They are httle sensitive to 
metallicity and to the SFH provided the SFR has been nearly constant in the last 10*^ years. 
Very young starburst galaxies would significantly deviate from a constant SFR model and 
require a higher value of rjuv, while the UV emission of early-type galaxies is contaminated by 
older stars. For an optically selected sample with a mix of SFHs and metallicities similar to 
the SDSS/GALEX sample defined in the previous section (the average metallicity of which is 
O.SZq), Samir et al. (2007) suggest using their median conversion factor \og(J]p^y) = 28.14. 
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We assume this calibration for both UV bands in the following. 

Figure [2] shows SFRuv against SFR^ for the NUV and FUV bands (left and right 
panel respectively), assuming no dust correction for the UV luminosities. The dashed green 
line denotes equality of SFR. Blue and red dots distinguish between blue and red galaxies 
defined as having ^■^{n — r) < 4 and > 4 respectively. The histograms in inset show the 
distribution of the SFRjjv to SFR^ ratios for the blue and red populations. As expected 
from uncorrected luminosities, SFRuv underestimates the 'true' SFR, including for the red 
population although part of their UV luminosity is unrelated to the current SFR. The scatter 
is large in both bands, indicating a large range of UV attenuations for a given SFR or a given 
UV luminosity. There is also a clear trend with SFR^ in the sense that g alaxies with higher 



SFR tend to require a larger dust correction, as was first observed by IWang fc Heckman 



fll996h 



3.2. Dust attenuation estimates 



The fraction of UV flux emitted by n ew born sta rs and absorbed by dust in the galaxy 
is reradiated at infrared (IR) wavelengths (lBuatlll992l ). Assuming a standard extinction law 
and that the dust is heated by intrinsically young stellar populations, the FUV attenuation 
can be approximated by: 

Afuv = 2.5 X \og{nIRX + 1) (2) 

where IRX is the ratio of the IR to UV luminosities (the so called infrared excess) and 
corrects for the fraction of IR luminosity heated by older stars and by light bluer than the 
FUV band (IMeurer et al.lll999l ). Other relations were derived that yield very similar results 
(e.g. Buat et al. 2005). UV reddening as measured by the slope [3 of the UV continuum 
(/a oc A^) or a UV color, correlates with IRX in starburst ga laxies, as expected from a 
foreground screen of dust (IWitt et al.l Il992l : ICalzetti et al.l Il994l ) . Thus [3 or UV colors are 
often used to estimate Apjjy- The IRX//3 correlation was recently revisited and corrected 
by several authors for more 'normal' star-forming galaxies using GALEX da ta (Seibert et 



al. 2 005, Cortese et al. 2006, Salim et al. 2007, Johnson et al. 2007). ISeibert et al. 



( 120051 ) (hereafter Se05) found that the starburst relation systematically overestimates the 
FUV attenuation of more quiet galaxy types by 0.58 mag, albeit with a large scatter. They 
derived the following empirical relation from a diverse sample of ~ 200 galaxies with UV 
photometry from GALEX and FIR photometry from the Infrared Astronomical Satellite 
{IRAS): 



A 



5e05 
FUV 



3.97 {mpuv - mNuv - 0.1) + 0.14 



(3) 
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where mpuv and rriNuv are the apparent magnitudes in the FUV and NUV bands respec- 
tively. The 0.1 magnitude offset corrects for a change in cahbration between the GALEX 
photometry used by Se05 (the internal data release IR0.2) and that used in the present paper 
(IRl.l) (Seibert et al., in preparation). 



Salim et al.l (120071 ) (hereafter Sa07) derived a yet shallower relation between the atten- 
uation and the UV color of normal blue galaxies using a different technique and the much 
larger GALEX/SDSS sample we are using here. They obtained the SFR and FUV atten- 
uations, among other physical properties, by fitting the UV and optical photometry to an 
extensive library of m odel SEDs for which dust attenuation was computed from the same 



Chariot fc Falll (120001 ) model used by B04. While B04 obtained their SFRs and attenuations 
from optical emission lines (the Ha line and the Balmer decrement to first approximation), 
Sa07's are essentially constrained by the UV fluxes and the UV colors respectively. The agree- 
ment is generally good between the two approaches but discrepancies remain, in particular 
between the attenuations as we'll discuss further below. Sa07 derived the following simple 
prescription from their extensive modeling for galaxies with °(n — r) < 4 and °(/ — n)< 0.9: 

.4|"t°^ = 2.99 °(/-n) + 0.27 (4) 

where the subscript refers to rest-frame colors fc-corrected to z = 0. The small (< 5%) 
fraction of galaxies with "(n — r) < 4 and "(/ — n)> 0.9 are assigned a constant attenuation 
of 2.96. 

Longer baseline colors such as UV— optical colors carry mixed but separable information 
about the SFH and IRX. Using a sample of galaxies with UV through IR photometry from 
GALEX, SDSS and Spitzer, Johnson et al. (2006, 2007) (hereafter J07) showed that given 
the SFH of a galaxy (they used D„(4000)), IRX could be more accurately inferred from 
UV— optical colors than from UV colors. Assuming Eq. [2] with /i = 0.6, they derived the 
following relation for galaxies with D„(4000)< 1.6: 

Af^y = 1.21 - 2.04x + 1.45?/ - 0.98?/^ (5) 

where x = D„(4000) - 1.25 and y = ° - r) - 2. 

We note SFRpuv,corr the SFR derived from the UV luminosity corrected for dust at- 
tenuation using one of the above equations. Figure [3] shows SFRFuv,corr against SFR^ for 
the blue galaxies (°-^(n — r) < 4) using Eq. [3] (Se05), H] (Sa07) and O (J07) as indicated. 
In each panel the dotted, dashed and solid blue lines show the ordinary least-square (OLS) 
regression of the Y axis on the X axis, t he OLS regressioii of the X axis on the Y axis and 



the bisector of those 2 lines respectively (llsobe et al.lll990l ). We choose the bisectors as the 



'best-fit' lines, here and in the rest of the paper. The best-fit slopes, variances, correlation 
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coefficients and residual scatters are listed in Table 1 for the 3 attenuation models (first 3 
lines). The histograms show the distributions of the SFRpuv,corr/ SFR^ ratios compared to 
the distribution of the uncorrected SFRpuy to SFR^ ratios. The averages of the distribu- 
tions are also listed in Table 1. All 3 methods provide a very good average correction with a 
reduced scatter compared to the uncorrected SFRpuv, especially so for the J07 correction. 
However a residual trend with SFRg remains in the sense that galaxies with the highest and 
lowest SFR tend to be under and over-corrected respectively. This indicates that the models 
do not quite sufficiently scale with the SFR to straighten up the uncorrected correlation in 
Fig. [21 The trend is minimal for the Se05 correction for which the scatter is largest and 
most pronounced for the J07 correction for which the scatter is otherwise best reduced. This 
trend with SFRg is the same as that noted by Sa07 as a trend with mass (their Fig. 8) and 
by J07 as a trend with '^■^{n — r) and D„(4000) (their Fig. 12 and 13). Indeed mass and to 
a lesser extent colors and D„(4000) correlate with SFR^.. The reason for it remains unclear 
but Sa07 concluded that the most likely interpretation in the framework of their modeling 
was that attenuations were less well constrained by the UV data than by the emission lines 
at the two ends of the distribution. In the case of J07, the parametric relation between IRX, 
D„(4000) and ^•^{n — r) is a good fit to blue galaxies, more so than between IRX and UV 
color (see J07 for a detailed discussion). Therefore it is perhaps the relation between IRX 
and Apuv that is not totally adequate. We return to this point in the next section. 



4. Reconciling UV and emission line Star Formation Rates 

4.1. Empirical color corrections 

Assuming as we do that SFR^ is the current best dust-corrected SFR estimate and that 
our choice of rjuv is adequate, the UV attenuation (FUV or NUV) can be directly measured 
as: 

Auv = -2.5 \og{SFRuvlSFR,) (6) 

We now revisit the color dependence of these known attenuations. Figure HI left panel, shows 
the SFRnuv to SFR^ ratios as a function of °'^(n — r) for the SDSS-fNUV sample. The 
solid lines show the OLS bisector for each D„(4000) bin as shown in inset. The 3 bins making 
up the blue sequence (D„(4000)< 1.6) add up to form a single tight correlation while the 
2 bins with D„(4000)> 1.6 form a scattered cloud. For galaxies with D„(4000)< 1.6 and 
0.1^^ — r) < 4, the bisector fit is: 



ANuv,n-r = 1-71 "'^ (n - r) - 2.86 



(7) 
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with a linear correlation coefficient r = 0.77 and rms = 0.70. The SDSS+NUV sample 
allows us to see the impact of the SFH on the attenuation/color relation. It is bimodal 
to first order: attenuation is linearly dependent on color for blue sequence galaxies and 
practically independent of it for red sequence galaxies. Similar correlations are found when 
using FUV luminosities and/or other UV— optical colors. Using the FUV luminosity and 
'■\f~g) yields: 

AFuv,f-a = lM'-\f-g)- 2.57. (8) 

The correlation is tighter than Eq. [7] (r = 0.84, rms = 0.66). On the other hand ^'^{f — n) 
results in a poorer and more scattered correlation (r = 0.56, rms = 1.02): 

Afuvj-u = 4.05 - n) - 0.18. (9) 

The last two correlations a re shown in the righ t panel of Fig. HI The latter correlation is 



similar to that proposed by ICortese et al.l (120061 ) (their Fig. 10) using log^L-na/ Lpuv) as a 



function of {mpuv ~^nuv) for a small sample of star-forming galaxies in the COMA cluster. 
It is consistent with Eq. [3](Se05) but a much steeper function of UV color than Eq. |l](Sa07). 
However all 3 equations as well as Eq. [S] converge for the majority of galaxies around the 
peak of the color distribution {^■^{n — r) ~ 2.4 — 2.5) and yield similar average values in good 
agreement with the measured average attenuation (< Apuv >~ 1-8). Equation [5] (J07) 
yields a slightly lower average attenuation of 1.6. Discrepancies between these corrections 
are largest for the reddest and bluest galaxies. 

We note SFRfuv,c the FUV based SFRs corrected using one of the functions of color 
derived in this section. Figure [5] shows SFRfuv,c against SFRe using Apuvj-n (left panel) 
and Apuv,f-g (right panel) for the blue population. The correlation is very close to equality 
in both cases, as expected since the corrections were designed to minimize SFRpuy^J SFR^, 
and the scatter is best reduced using '^'^{f—g), as expected as well from the higher correlation 
coefficient in Eq. [HI However a small residual trend with SFRe remains in this case, which 
means that the SFR dependence of the attenuation is not completely accounted for by the 
color dependence {\og{S F Ruv / S F R^) would be better fitted by a linear function of color 
plus a linear function of \og{SFRg)). The trend practically disappears when using °'^(/ — n) 
but the correlation is less significant, as in the case of Se05. The parameters of the fits are 
listed in Table 1 for comparison with the corrections presented in the previous section. There 
is a small trade-off between the scatter and the trend with SFR^ (the lower the rms, the more 
the slope deviates from unity) except for the empirical °'^(/ — g) correction for which the 
combination of trend and scatter is best reduced (a > 0.9 and rms < 0.3). Figure E] shows 
the difference between the measured attenuation (Eq. [6]) and the 4 parametric estimates 
(Se05, Sa07, J07 and Eq. [8]) as a function of SFRe {AApuv = Apuv — ^FC/v,modei)- The red 
curves are isodensity contours. The dashed green lines mark AApuy = — 1, and 1. All 4 
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methods converge with Apuv around the peak of the SFR distribution (~ 2M0yr~^ ) and as 
noted above provide good average corrections but Eq. [8] minimizes AApuv for the majority 
of the galaxies as well as the dependence with SFR^. Although the improvement is by no 
means dramatic, it provides an estimate of the FUV attenuation which best recovers the 
SFR derived from emission lines at an equivalent or lower 'cost', since at least one optical 
photometric measurement and the redshift are required for fc-corrections in all cases. 



4.2. Discussion 

As aperture corrections are an important source of uncertainty in deriving total SFRs 
from the SDSS fiber spectra, we check that the above correlations are not affected by aperture 
effects. We define the aperture correction (AC) as the ratio of the total SFR (SFRe) to the 
SFR estimated within the fiber (B04) and split the blue sample into 3 bins of increasing 
aperture correction. Figure [7| shows the SFRpuv to SFR^ ratios as a function of °'^(/ — g) 
for the 3 bins (shown in inset), shifted by a constant as indicated for clarity. The solid 
lines show the bisector fits in each bin and the dashed lines show Eq. [8] shifted by the 
appropriate amount for comparison. No significant difference is seen in the correlation itself 
as a function of aperture correction, but galaxies with SFRpuy > SFR^. (yielding a negative 
FUV attenuation) have systematically high aperture corrections (and low redshift). It is 
likely that in these large nearby galaxies, the fiber missed off-centered regions of enhanced 
star- formation and that SFR^ was underestimated. 

Another source of uncertainty may lie in the definition of the UV attenuation (Eq. [6]). 
Allowing Tjuv to vary, the SFRjjy to SFR^ ratio would no longer be a measure of attenuation 
alone but a combination of Apuy and 'quv'- ^og{S F Ruv / S F Rf,) = —OAAuv + ^og^rjuv/Vuv)- 
B04 showed that the conversion factor t^ho from Ha luminosity to SFR {rjna = L^a/ SFR^) 
decreases with mass (and metallicity), spanning nearly 0.4 dex in the range of mass (and 
metallicity) spanned by the SDSS sample. This is interpreted as massive/metal rich galaxies 
producing less Ha than low mass, metal poor galaxies for the same SFR. However rjuv is 
much less sensitive to metallicity than Ha, consequently SFRuy/SFRe is expected to be 
a nearly direct measure of UV attenuation except perhaps for galaxies with very low UV 
attenuation and metallicity. In any case the ~ 1 dex variation in rjuy required to straighten 
up the trend seen in Fig. El most notably for the JOT dust correction, is definitely ruled 
out. Let's note that Eq. [Hlcan be used to recover SFR^ from the observed UV luminosity 
whatever the interpretation of the SFRuy/SFRe ratio provided rjuy is assumed in Eq. [1] 
Figure [8] shows the relation between Ajquy and Apuv as defined in Eq. [61 The solid blue 
line is the best linear fit; the dashed green line at Ajsiuv = 0.75 Apuv is the ratio expected 
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from a A absorption curve (ICharlot &: Fal]||2000l ). The distribution of the A^uv to Apuv 



ratios is shown in inset as a sohd histogram, and that of the A^uy^n-r to Apuvj-g ratios as a 
dotted histogram. The median ratio for the measured quantities is 0.74 (0.75 for the fits), in 
excellent agreement with the prediction. The GALEX data are therefore consistent with the 
attenuated UV fluxes predicted by the dust model used to derive SFRe from independent 
emission line measurements. 

As noted in the previous section, the relation between IRX and Apjjy may be responsible 
for the discrepancy between Apuy and AjP^y (Eq. \5^. Figure [HI shows the relation between 
Afuv (Eq. ED and IRX derived from D„(4000) and °-^(ra - r) (Eq. [2] and [5]). The sohd line 
is the relation used by J07 to relate the two quantities (Eq. [2] with = 0.6). While it is 
appropriate for the average galaxy, it becomes discrepant at the blue and red ends (or for the 
least and most massive/star-forming galaxies). There is also a large range of attenuations 
for a given IRX. A color (or SFR) dependent fi parameter accounting for the fraction of IR 
flux due to new born stars (as opposed to preexisting older stars) would remedy some of the 
discrepancy. A blue galaxy would have little dust attenuation, little FIR emission for a given 
UV flux (a small IRX) and its little FIR emission would seem to have little to do with the new 
born stars (small /i). At the other end, a red (star- forming) galaxy would be a dusty galaxy 
with a high SFR and its large FIR emission (large IRX) would be entirely due to the heating 
of dust by its new stars (large fi). A fi parameter as large as 2 (meaning that the obscured 
UV emission would have to be twice the observed FIR emission) is actually necessary to 
reach the upper enveloppe of the IRX/ Apuy distribution. This might be accounted for by 
the uncertainty in the IR flux estimate, or by a large fraction of obscured UV photons being 
reprocessed at wavelengths other than IR. J07 modeled the contribution of new born stars 
to IRX as a function of galaxy color and arrived at the opposite conclusion: red galaxies 
have a higher contribution of older stars to their IR emission than blue galaxies and should 
therefore required a lower n, making the trend between the corrected UV luminosity and 
SFRg even more pronounced. The validity of SFRe at low and high mass may of course 
be questionned but it seems that the estimate of the IR flux and the interpretation of IRX 
in terms of FUV attenuation currently involve more uncertainties than the interpretation of 
the optical data. 



4.3. Limitations: the oldest and youngest galaxies 

The color corrections do not apply to red sequence galaxies in the local sample. Dust 
attenuation estimates based on correlations between IRX and colors do not apply well to 
early-type galaxies either, both the UV and IR SEDs of such galaxies being much less directly 
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related to the emission of young stars than those of late- type galaxies (J07). This isn't a 
drastic limitation to the various methods since red sequence galaxies contribute little to the 
overall SFR in the local Universe, and even less as redshift increases. But we would like 
to know whether the above correlations between attenuation and UV-optical color apply to 
star- forming galaxies at higher redshift, where rest-frar ne UV fluxes are ge nerally corrected 
using the IRX//3 correlation of local starburst galaxies (IMeurer et al.lll999l ). 



Although no spectroscopic data exist at high redshift that allow the same emis s ion line 
fitting technique as the SDSS spectra, we can use the unique sample of lErb et al.l (l2006bl ) 
who were able to acquire Ha flux measurements for 114 UV selected galaxies at z ^ 2, for 
which optical and NIR photometry is also available. They defined SFRfjjy as in Eq. [1] 
a nd SFR^ rv = -^H r>/?7Hn-, using the rjpuv and ?7Ha values of iKennicuttI (119981 ) converted to 



Chabrierl ( 20031 ) IMF q A factor of two aperture correction was also applied to the Ha 



luminosities. Dust corrections were derived from the best-fit values of E{B — V) obtained 
from fitting SEP ni odels to the mu l tiband photometry and using the extinction law of 
Calzetti et al.l (120001 ) (see lErb et al.l (l2006bl ) for details). The authors assumed that the 
color excess of the nebular emission lines was equal to t hat of the UV continuum, rather 
than 2.5 times larger as proposed by lCalzetti et al.l (120001 ) . as it yielded the best agreement 
between the UV and Ha SFRs after dust correction. The corrected SFR^a are at most 3 
times the uncorrected values and less than twice for most of the sample. 

As high redshift galaxies have lower metallicities than local galaxies on average, a higher 
value of 7]iia might be justified for this sample. However as both the range of masses and the 
range of metallicities span by the z = 2 galaxies remain within th ose of the SDSS sample 
despite evolution in the mass/metallicity relation (lErb et al.ll2006al ). rjua is not expected to 
be larger than the value predicted for the least massive/most metal poor galaxies in the 
local sample, i.e. a factor of 1.5 higher than the Kennicutt value used by Erb et al. (B04, 
their Fig. 7). We assume the uncorrected values of SFRna divided by 1.5 to be lower limits 
to the SFR and use the dust corrected values as upper limits. We computed the absolute 
magnitudes of the galaxies in the GALEX and SPSS ban ds from their optical and NIR 
photometry using kcorrect v4_l ( IBlanton fc Roweid 120071 ) . 

Figure [TUl shows SFRj^uy / SFR^. as a function of ^-^{n — r) assuming SFRna = SFR^ 
for the high redshift sample and rji^uv = Vfuv all galaxies as before. The pink circles 
are the uncorrected values of SFR^a- The upper and lower limits are defined as above. The 



^To do so, Erb et al. multiplied the Kennicutt factors (computed for a Salpeter IMF) by 1.8, however 
a conversion factor of 1.58 between the Salpeter and Chabrier IMFs is more appropriate (S. Salim, private 
communication). We corrected their SFRs accordingly and multiplied them by 0.94 to account for the small 
difference between the Chabrier and Kroupa IMFs. 
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local blue population is plotted in blue with the correlation derived in the previous section 
(Eq. [7]). It is clear that the majority of LBGs which have extremely blue colors do not follow 
the same correlation as the local galaxies, but the reddest ones ("•^(ra — r) > 1.5) may still 
be consistent with the local correlation or show a similar trend with color with only a small 
blue shift. The bluest LBGs cluster blueward of the local correlation in a region of low UV 
attenuation. 

Also overplotted is a sample of 97 compact UV Luminous Galaxies (UVLGs) drawn 
from the present sample and from a cross- match betw een the SDSS and the larger, shallower 
GALEX All Sky I maging Survey 



of them. UVLGs (IHeckman et al. 



2005 



Hoopes et al. 2006). SFR^ from B04 are available for all 



Hoopes et al.ll2006l ) are locally rare galaxies defined 



as having FUV luminosities typical of LBGs: Lpuv > 2 x 10 L(7^, corres ponding to ~ 0.3L 



at z ~ 3 fISteidel et al.lll999f ) but to ~ 5L^ at z ~ flWvder et al.l 120051 ). While low surface 
brightness UVLGs are simply extra large versions of normal spiral galaxies, high surface 
brightness UVLGs with Ipuv > lO'^L0kpc~^ were found to consist primarily of compact 
starburst systems. Among these, the "supercompact" UVLGs with Ipuv > lO^Lokpc"^ 
bear a remarkable resemblance to high redshift LBGs for a wide range of physical properties 
(mass, SFR, rnetalli city). They are thought to be their closest analogs in the Local Universe 
( iHoopes et al.l 120061 ). Compact and supercompact UVLGs are represented in Fig. [TO] with 
open squares and filled triangles respectively. Both categories occupy the same region of the 
plot as the high redshift sample. The supercompact UVLGs are unusually blue {^■^{n — r) < 
1.5) among the local galaxies and like the bluest LBGs, lie the furthest away from the 
bulk of the local population. The extreme blue colo rs are an indicatio n of strong recent star- 
formation (as well as low attenuation). As noted by iKennicuttI (119981 ). the calibration rjuv in 
Eq. [1] might be significantly higher for strong starburst galaxies such as these. A higher rjuv 
would lower their SFRuy to SFR^ ratio proportionally and further separate them from the 
main population. Therefore very blue galaxies with very recent star-formation, both locally 
and at high redshift, form a distinct cluster of their own blueward of the attenuation/color 
relation of the blue sequence. This locus adds to the blue shift of the attenuation/color 
relation with D„(4000) seen in Fig. H] (left panel) as those galaxies would have D„(4000) 
indices lower than our lowest bin. The reddest of the compact UVLGs and of the LBGs 
cover a wider range of attenuations which seem to correlate with colors in the same way as 
local galaxies only shifted to bluer colors. 

D„(4000) estimates are not available for the high redshift sample and only available for 
a fraction of the UVLGs but we can use the inverse of the specific SFR - Tsfr = M^^/SFR 
- as a common SFR time-scale for the lo cal and high reds hift samples. Stellar masses were 
estimated for most of the z ^ 2 saniple (|Erb et al.ll2006d ). Masses for the GALEX/SDSS 



sample are from iKauffmann et al.l (l2003al ). A third of the UVLGs have mass estimates 
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from this catalog; for th e remaining 2 /3 we use the values derived by iHoopes et al.l (120061 ) 
via SEP fitting following ISalim et al.l (I20051) . Th e agreement between the mass estimates of 
Kauffmann et al.l (l2003al ) and lSalim et al.l (120051 ) for the GALEX/SDDS sample is good with 
a rms of 0.12. As above, we use SFR = SFR^ for the local sample including the UVLGs and 
SFRna for the 2 ~ 2 galaxies. Figure [TT] shows Tsfr as a function of °-^(n — r). The SDSS 
galaxies are color coded in bins of D„(4000) as shown in the inset of Fig. H] (this highlights 
the relation between D„(4000) and Tsfr)- The LBGs are represented by filled circles, the 
compact UVLGs by open squares and the supercompact UVGLs by filled triangles. The 
two horizontal lines correspond to the age of the Universe aX z = and 2 (13.8 Gyrs and 
3 respectively). Roughly, galaxies with Tsfr larger than the age of the Universe at their 
redshift (e.g. nearly all local galaxies with D„(4000)> 1.6 and a few LBGs) have had larger 
SFR in the past. Inversely, galaxies with Tsfr shorter than the age of the Universe at 
their redshift (most LBGs and nearly all the compact UVLGs) must be forming stars more 
intensely than in the past. 

Figure[TT]is an analog of Fig. [T]with the addition of an 'ultrablue' sequence at °-^(n— r) < 
1.5 and Tsfr < 3 Gyrs consisting of young compact starburst galaxies. The Auv /color 
correlations derived in the previous section hold for galaxies with rather uneventful SFHs (3 
Gyrs ^TsFR ^ 15 Gyrs). They are a majority today but may not be when the Universe was 
only ~ 3 Gyrs old, although many LBGs look very much like local blue sequence galaxies. 
Furthermore a dominant fraction of the stel lar mass at z > 2 is found in redder galaxies 



which are largely abse n t from UV s urveys (iRudnick et al.l l2006t Ivan Dokkum et al.l 12006 



Marchesini et al.ll2007l ). iKriek et al.l (l2006bl ) showed that almost half of their sample of NIR 
selected galaxies at z = 2.0 — 2.7 have low SFRs and Tsfr > 10 Gyr (from their Fig. 2). 
These galaxies are redder than LBGs and would lie in the same part of the plot as the local 
population. 

Although we can't conclude on the use of the local Auy / color correlation at high redshift, 
we may expect it to hold for most galaxies to at least intermediate redshifts or to be slightly 
shifted to the left as only 'mild evolution' of the blue sequence was reported between z = 
and 1 (no change in the number density and colors only ~ 0.3 mag bluer; Blanton 2006 ). The 
luminosity density from UVLGs was found to undergo dramatic e volution between z = an d 
1, reaching > 25% of the total FUV luminosity density at z = 1 (ISchiminovich et al.ll2005l ). 
but this includes all UVLGs, i.e. mostly very large but otherwise ordinary spiral galaxies. 
Compact and supercompact UVLGs are very rare in the local Universe. Although their 
evolution with redshift is yet unknown, they are unlikely to dominate the galaxy population 
at z = 1. 
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Conclusions 



Using a large sample of galaxies from the Sloan Digital Sky Survey spectroscopic catalog 
with measured SFRs and UV photometry from the GALEX Medium Imaging Survey, we 
derived empirical linear correlations between the UV attenuation measured by the SFR to 
observed UV luminosity ratio, and the UV— optical colors of blue sequence galaxies 



(0.1, 



n — 



r) < 4). The SFRs were derived from a detailed modelling of the emission lines in the optical 
spectra (IBrinchmann et al.ll2004j ) and were considered best estimates. The attenuation/color 
relation provides a simple prescription to correct UV measurements for dust attenuation 
in the absence of SDSS quality data. We found or confirmed that other UV attenuation 
estimates (Seibert et al. 2005, Salim et al. 2007, Johnson et al. 2007) tend to over (under) 
correct the UV luminosit y of galaxies with the lowest (highest) emission line SFRs or mass. 
Using a sample of LBGs (lErb et al.ll2006bl ) at z ^ 2 with measured Ha emission as well as 



a sam ple of local compact UV luminous galaxies with LBG like properties (IHoopes et al. 



20061 ) . we found that extremely blue galaxies at both low and high redshift escaped the 
attenuation/color relation of the blue sequence to form a low attenuation sequence of their 
own. As such g alaxies are very rare locally and th e blue sequence does not evolve much 
from z = to 1 (jBlantoru l2006l : IWillmer et al.ll2006l ). we expect our attenuation correction 
to remain adequate for the majority of galaxies to at least intermediate redshifts. 
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Table 1. Correlation parameters between the corrected FUV SFR {SFRpuv,modei) and 
SFRe for the various attenuation models described in the text as hsted in the first column. 
The next columns list the slope a of the linear correlation fitted to the SFRFuv,modei vs 
SFRe relation, the variance of the slope, the correlation coefficient, the rms and the 
average of the SFRFuv,modeil SFRe ratios. 



Model 


slope a 


Var(a) 


r 


rms 


< SFRpuV-model/ SFRe > 


Se05 


0.9498 


5.181C-05 


0.7392 


0.4209 


0.0012 


Sa07 


0.8314 


2.23e-05 


0.7866 


0.3328 


0.0421 


J07 


0.7171 


1.242e-05 


0.8658 


0.2278 


-0.0847 


Eq. 9 


0.9799 


5.14e-05 


0.7487 


0.4266 


0.0061 


Eq. 8 


0.9258 


1.915e-05 


0.8705 


0.2893 


0.0316 
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Fig. 1.— The 4000A break index, D„(4000), as a function of °-i(n - r) for the SDSS+NUV 
sample. The blue curves represent lines of equal FUV attenuation by Johnson et al. (2007). 
The horizontal and vertical dotted lines at D„(4000)= 1.6 and ^■^{n — r) = 4 mark the 
boundaries between the 'red sequence' and 'blue cloud' populations. 



- 23 - 




Fig. 2. — The uncorrected U V SFRs (SFRmrv o'^^ th e left; SFRpjjy on the right) against the 
emission hne derived SFRg (iBrinchmann et al.ll2004l ) in units of Mq yr^^. The dashed green 
hues shows SFR equahty. The blue and red dots represent blue and red sequence galaxies, 

n — r) > 4 respectively. The blue and red histograms in 



defined as ^'^(n — r) < 4 and 



0.1, 



each panel are the \og{SFRuv/ SFRg) distributions for the two populations. As expected 
from uncorrected luminosities SFRuv underestimates SFRe, especially so at high SFR. 
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Fig. 3. — The dust corrected FUV SFRs {S F Rpuv,corr) again st SFRp for t he blu e popu- 
lation 



(O.l, 



n 



)< 4) using th e FUV attenuation m ethods of ISeibert et al.l (120051 ) (Se05), 
Salim et al.l (120071 ) (Sa07) and IJohnson et al.l (120071 ) (J07) as indicated. The SFRs are in 
units of Mq yr~^. The dotted, dashed and solid blue lines in each panel are the ordinary 
least-square (OLS) regression of the Y axis on the X axis, the OLS regression of the X 
axis on the Y axis and the bisector of those 2 lines respectively. The histograms show the 
distributions of the SFRpuVfiorrl SFR^ logarithmic ratios: light blue, red and orange for 
Se05, Sa07 and J07 respectively. The dotted histogram is the distribution of the uncorrected 
SFRpuY / SFRe ratios. All three models provide a very good average correction but tend to 
over (under) correct galaxies with the lowest (highest) SFRs. 



- 25 - 




Fig. 4. — The color dependence of the SFRuv to SFR^ ratio for different bins of 
D„(4000). Left: \og{SFRNuv / SFRe) against - r) for the SDSS+NUV sample; Right: 
hglsFRpuv/SFR^) against - g) and against - n) for the SDSS+NUV+FUV 
sample. The data are color coded to match the binned D„(4000) distribution histogram 
shown in inset. The colored solid lines show the fitted linear correlations in each bin of 
D„(4000). The 3 bins making up the 'blue cloud' population (D„(4000)< 1.6: light blue, 
dark blue and green) are tightly correlated with color, while the 2 bins forming the red 
sequence (D„(4000)> 1.6: yellow and red) are practically independent of it. The color cor- 
relations observed for blue galaxies provide simple UV attenuation corrections (see text for 
details). 
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Fig. 5. — The color corrected FUV SFRs against SFRf,. The FUV luminosities are corrected 
using the correlations established in Fig. H] between Apjjy and °'^(/ — n) (Eq. [U left panel) 
and between Apuv and ^'^{f — g) (Eq. [HI right panel). The histograms in inset show the 
distributions of log^SFRpuv/ SFRe) with and without corrections (solid and dotted hues 
respectively). SFRe is best recovered using the Apuv/^'^if ~ d) correlation. 
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Fig. 6. — AApuv = AFuv—AFUv,raodei as a function of SFRe for the various models described 
in the text as indicated (Sa05, Se07, J07 and Eq. [S]). The dashed hues mark AApuv = 
and +1. The red curves are isodensity contours. All 4 methods converge with Apjjy around 
the peak of the SFR distribution (~ 2M0yr^^ ) and provide good average corrections but 
Eq. E] (the Apuv/^-^if - g) correlation) minimizes AApuv and its dependence with SFR^ 
for the majority of the galaxies. 
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Fig. 7. — The SFRpuv to SFR^ ratio as a function of °'^(/ — g) va. ?> bins of aperture 
corrections (AC) as shown in inset. The AC is defined as the ratio of the total SFR {SFRg) 
to the SFR within the SDSS fiber. For clarity the data are shifted by 0, +1 and +2 for the 
first, second and third bin respectively. The solid lines are the linear fits in each bin while 
the dashed lines show the correlation for the full sample (Eq. [8]) shifted by the appropriate 
amount. Aperture effects do not bias the correlation. 
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Fig. 8. — The relation between A^uv and Apuy as defined by the SFRuy / SFRe ratios (Eq. 
|6]). The sohd blue fine is the fitted correlation; the dashed green line at Ajsiuv = 0.75 Apuv 
is the ratio expected from a A"'''^ absorption curve (Chariot & Fall 2000). The distribution 
of the A]^uy to ApuY ratios is shown in inset (solid histogram) with that of the fits ratios 
{A]^uv^n-r/ Apuvj-g] dotted histogram). The median ratio for the measured quantities is 
0.74 (0.75 for the fits), in excellent agreement with the prediction. 
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log(IRX) 



Fig. 9. — The FUV attenuation (Eg. El) as a function of IRX derived from D„(4000) and 
0.1 — r) following Johnson et al.l (2007h for the blue population. The solid line is the 
^FC/v/IRX relation of iMeurer et al.l (jl999) (Eq. [2]) assuming n = 0.6. The lower and upper 
dotted lines are for /x = 0.1 and 2 respectively. 
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Fig. 10. — The SFR^uv to SFR^ ratio as a function of — r) for the local blue popu- 
lation (blue dots), the sample of local compact and "supercompact" UV Luminous Galaxies 
(UVLGs) of Hoopes et al. (2006) (open square s and filled yellow triangles respectively) and 
the 2; ~ 2 Lyman break galaxy (LBG) sample of lErb et al.l (j2006bl ) assuming SFR^a = SFR^ 
(pink filled circles; see text for details). The straight line is the correlation fitted to the local 
blue population (Eq. [7]). 
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Fig. 11. — The SFR timescale Tsfr = Mj,/SFRe as a function of °'^(n — r). Local galaxies 
are shown by dots color coded in bins of D„(4000) as in Fig. |H UVLGs and LBGs are 
represented as in Fig. [TOl The two horizontal dashed lines mark the age of the Universe at 
z = 2 and 2; = (3 and 13.8 Gyrs respectively). This plot is similar to Fig. [1] with the 
addition of an ultrablue sequence at ^-^{n — r) < 1.5 and Tsfr < 3 Gyrs consisting of young 
compact starburst galaxies for which the local UV attenuation/UV-optical color relations 
derived in this paper do not apply. 



